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ABSTRACT 

We analyze intermediate degree p-mode eigenfrequencies measured by GONG 
and MDI/SOHO over a solar cycle to study the source of their variability. We 
carry out a correlation analysis between the change in frequencies and several 
measures of the Sun's magnetic activity that are sensitive to changes at different 
levels in the solar atmosphere. The observations span a period of about 12 
years starting from mid- 1996 (the minimum of cycle 23) to early-2008 (near 
minimum of cycle 24), corresponding to a nearly complete solar activity cycle. We 
demonstrate that the frequencies do vary in phase with the solar activity indices, 
however the degree of correlation differs from phase to phase of the cycle. During 
the rising and declining phases, the mode frequency shifts are strongly correlated 
with the activity proxies whereas during the high-activity period, the shifts have 
significantly lower correlation with all activity proxies, except for the 10.7-cm 
radio flux. In particular, the proxies that are only influenced by the variation 
of the strong component of the magnetic field in the photosphere have a much 
lower correlation at the high-activity period. On the other hand, the shifts are 
better correlated with the proxies sensitive to changes in the weak component of 
the magnetic field. Our correlation analysis suggests that more than 90% of the 
variation in the oscillation frequencies in all activity phases can be explained by 
changes in both components of the magnetic field. Further, the slopes obtained 
from the linear regression analysis also differ from phase to phase and show a 
strong correlation with the correlation coefficients between frequency shifts and 
solar activity. 

Subject headings: Sun: helioseismology - Sun: oscillations - Sun: activity - Sun: 
magnetic fields 
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Introduction 



Understanding the source of the varying solar oscillation eigenfrequencies has been a 
topic of interest since the mid-eighties. A decrease in low-degree mode frequencies be- 
tween 1980 (near so l ar ma ximum) and 1984 (near solar minimum) was init ially reported by 



Woodard fc Noyes I (119851 ) . Later, using Big Bear Solar Observatory data, IWoodard et al 



( 119911 ) showed that the frequencies were strongly correlated with the magnetic field strength. 
The variation of the frequencies with the changing magnetic activity is now well estab- 



lished. Suc h changes affect bo t 



l ow-degree (ISalabert et al 



2004: 



i the central frequ encies and the frequency splittings of 



Chaplin et a. 



120071) as well as intermediate-degree modes 



feachmann fc Brownlll 9*931 : 1 Jain fc Bhatnagaii2003l : lAntiall20~03l : bziembowski fc Goodehoosl ) 
and exhibit a good correlation between the mode frequencies and the solar activity. The 
high-degree mode frequencies calculated using the Dynamics Program data of the Michelson 
Doppler Imager (MDI) on b oard Solar and Heliospheric Observatory fSOHO) also clearly 
show a solar cycle variation ( Rabello-Soares. Korzennik. &: Schou 20061 ). 



Detailed studies have shown that the mode frequencies do not al ways follow the same 



paths during the rising and the declining phases of the solar activity (jJimenez-Reyes et al. 



19981 : iTripathy et al.l 120011 ). In addition, using the continuous oscillation da ta from the 



Globa l Oscillation Network Group (GONG) for the rising phase of cycle 23, iHowe et al. 



( 119991 ) demonstrated that, although the average frequency shift is correlated well with mea- 
sures of surface magnetic activity integrated over the whole disk, the even-order splitting 
coefficients had a higher correlation with the corresponding coefficients of a Legendre polyno- 
mial decomposition of the surface magnetic field than with the total flux. A recent analysis of 
the year-wise distribution of the frequency shifts with the change in activity indices showed 
that bot h the linear-regressio n slopes and the magnitude of the correlation varied from year 
to year (ITripathy et al.l 120071 ). Their analysis further demonstrated that there were signifi- 
cant differences between long-term and short-term variations and concluded that averages of 
short-term variations do not ac curately reflect the long-term variations and vice versa. The 
analysis of I Chaplin et al.l (120071 ) for low degree modes involving three solar cycles also indi- 
cate that the frequency shifts are primarily sensitive to the weak component of the magnetic 
flux. 

Thus, the relationship between oscillation frequencies and solar activity indicators seems 
to be more complicated than the linear relation generally demonstrated in earlier studies 
( IWoodard et al. I Il99ll : iBachmann fc Brownl Il993l : iJain fc Bhatnagarl 120031 ) . With access to 
the high quality uniform data from GONG and MDI instruments for a nearly complete solar 
cycle, it is timely to revisit the temporal variability of the frequencies with the progression 
of solar activity cycle to understand the detailed relationship between them. In this paper, 
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we closely examine the relationship between the frequency shifts of the intermediate degree 
modes and nine activity proxies for cycle 23. We further divide the span of the activity cycle 
into three phases and compare the results obtained for different phases with the entire activity 
cycle. The frequency and activity data used here are presented in Section 2. We present 
results for individual data sets and their comparison in Section 3. A possible scenario to 
explain the variation in oscillation frequencies is discussed in Section 4 and the main findings 
are summarized in Section 5. 



2. Data 

2.1. Mode Frequencies 

The analysis presented here uses intermediate-degree mode data sets obtained from the 
GONG0 and MDJj, covering a period of about 12 years i.e. a nearly complete solar cycle 
starting from 1996 May 1. Each data set consists of centroid frequencies u Uj £ and splitting 
coefficients eij and each n, i multiplet is represented by a polynomial expansion 

Vn,t,m = Vn,£ + ^ Oifa, £)P>^(m), (1) 
i=l 

where P^ (m)'s are orthogonal p olynomials of degree i and i max is the number of a coefficients 



used in determining frequencies ( iSchou I Il992t ) .The remaining symbols in equation (1) have 



their usual meanings. The fitting is done using i max = 9 for GONG and 18 for MDI. 

For the study reported here, we cover a period from the minimum of solar cycle 23 (1996 
May) to the near minimum period of solar cycle 24 (2008 February). The continuous GONG 
data sets are centered at GONG Month (GM) 12 to 124 while MDI data have two gaps in 
1998-99 due to the breakdown of the satellite. It should be also noted that each 108-day 
GONG data set starts 36 days later than the previous one having an overlap of 72 days 
while the 72-day MDI data sets are non-overlapping. The MDI data have been obtained 
from the medium-^ program that provides almost continuous coverage to modes with i < 
300. Various other details of the data sets are given in Table I and duty cycles are plotted 
in Figure 1. We further note that the GONG data sets have relatively low duty cycle as 
compared to the MDI. The number of modes present in each GONG data set varies between 
1700-1800, but only one-third of these modes are observed in all 113 overlapping data sets. 



1 ftp://gong. nso.edu/data/ 

2 http:/ /quake. stanford.edu/ schou/anavw72z/ 
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Similarly, each MDI data set has approximately 1900 modes, which reduces to 885 common 
modes in all 58 data sets. The common mode sets in GONG and MDI are shown in Figure 
2(a-b). It should be noted that all these common modes are not successfully fitted for both 
projects, this might be due to the consequence of finite lifetime and stochastic nature of the 
modes. The common modes observed in both data sets are shown in Figure 2c. 



2.2. Proxies for Solar Activity 

To study the changes in mode frequencies with solar activity, we have used nine proxies 
of solar activity probing the solar atmosphere at different heights. Each activity index is 
briefly summarized below; 



Ri, the international sunspot number computed from the total number of spots visible 
on the face of the sun (JV) and the number of groups (G) into which they cluster. Rj 
is equal to 10G + N. 

KPMI, the line of sight magnetic field strength observed with 1 arc-sec pixels using 
the phot ospheric Fe I line at th e Kitt Peak Vacuum Telescope and averaged over the 
full disk ( iLivingston et al.lll976l ). 



MWSI, the Mt. Wilson sunspot index calculated using Mt. Wilson magnetograms 
taken at the 150-Foot Solar Tower. MWSI is determined by summing the absolute 
values of all pixels with a magnetic field strength (> 100 G) which is then divided by 
the total num ber of pixels (regardless of magnetic field strength) in the magnetogram 
Juirichlll99lh . 



F w , the integrated emission from the solar disc at 2.8 GHz i.e. 10.7 cm wavelength. It 
mainly represents the contributions from sunspots and radio plages in the upper chro 



mosp here in addition to the quiet-Sun background emission (jCovingtonlll969l ; iKundu 



19651 ). 



5. He I, the equivalent width of the He I 1083 nm solar absorption line. It is affected by 
photo-ionization radiation from the upper transition region and coron a. The index is 
calculated by measuring the change in temperature in transition region (jHarveyll 19841 ). 



6. MPSI, the magnetic plage strength index calculated in the same manner as the MWSI 
but the absolute value of the magnetic field strength is summed for fields between 10 G 
and 100 G. 
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7. Mg II, the Mg II core-to-wing ratio derived by taking the ratio of the h and k lines 
of the solar Mg II feature at 280 nm to the background or wings at approximately 
278 nm and 282 nm. The h and k lines are variable chromospheric emissions while the 
background emissions are more stable. This ratio i s a robust measure o f chromospheric 
activity, mainly for solar UV and EUV emissions (IViereck et al.ll200ll ). 



8. TSI, the total solar irradiance, describes the radiant energy emitted by the Sun over all 
wavelengths that falls each second on 1 square meter outside the earth's atmosphere. It 
is affected by the changes in both photospheric and chromospheric magn etic structures. 
We u se observations made by VIRGO instrument on board SOHO (IFrohlich et al. 



19970 . 



9. CI, the coronal index measures the total energy emitted by the sun's outermos t atmo- 
spheric layer (the corona) at a wavelength of 530.3 nm (IRybansky et al.lll994l ). 



It is worth mentioning that some activity indices are direct measurements of weak, strong 
or both components of the magnetic field at different heights in the solar atmosphere while 
others are surrogates of magnetic measurements. The strongest fields are concentrated in 
flux tubes in sunspots while the larger area around sunspots harbors active regions and their 
plages contain small-scale pores with relatively weak field strength. The weak components 
are more widely distributed in latitude than the strong components which are confined to the 
activity belts. Since flux tubes expand like a canopy in chromosphere above the photospheric 
footpoints, the field strength decreases with increasing height. 

Table 2 summarizes each activity index, including the spectral line of observation, the 
period of observation used in the analysis, the dominant component of the magnetic field, 
and the source of the data. A mean activity, /, is computed for each activity index over 
the same temporal intervals corresponding to the individual frequency measurements. This 
mean value is used in subsequent calculations. 



3. Analysis and Results 
3.1. Temporal variation in frequencies 

The centroid frequency shift 5v is calculated from the relation 

M*) = £%K<(*)/£%. (2) 
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where 8p n ^(t) is the change in the measured frequency with respect to a reference frequency, 
a Hj £ is the uncertainty in the fre quency measurement for a given n, £ multip let and Q n ^ 
is the inertia ratio as defined by IChristensen-Dalsgaard fc Berthomieu I (ll99ll ). There are 
several definitions for selecting a reference frequency, which have been discussed in detail by 
Howe et al.l (120021 ) . For a given mode, we adopt an approach of taking a temporal mean over 
all the available data sets. It has also been shown earlier that the frequency shift is strongly 
dependent on th e frequency and increases with increasing frequency (e.g. see Figure 3 of 
Jain et al.l (120001 ) ). We therefore restrict our analysis to common modes observed in all data 
sets to obtain an unbiased variation. 

We use 235 modes observed in all GONG and MDI data sets in the degree range of 21 
< i < 147, radial order range of 1 < n < 16 and the frequency range of 1.5 < v < 4.0 mHz. 
These modes are shown in Figure 2c. The reference frequencies are calculated by taking an 
average of the frequencies of a particular multiplet (n,£) from all of the available GONG 
and MDI data sets. Figure 3 shows the temporal variation of frequency shift for both data 
sets. The dashed line depicts the variation of the solar activity represented by radio flux, 
Fiq. It is clearly seen that the frequenci e s vary in phase with t he so lar activity as shown 
in earlier studies (e.g. IJain fc Bhatnagarl (120031 ); IChaplin et al.l (120071 )). We also note that 
there is a systemati c offset between G ONG and MDI frequency shifts, which is consistent 
with the findings of iHowe et al.l (120021 ) . The MDI frequency shifts are consistently higher 
than the GONG and this deviation is relatively large at the low activity periods. However, 
the overall shifts from the minimum to the maximum in both cases are approximately same. 
The discrepancy observed might be due to the fact that the GONG and the MDI pipelines 
consider time series of different lengths and use different mathematical methods to determine 
the mode frequencies. A detailed correlation analysis between the frequency shifts and solar 
activity is presented in the next section. 



3.2. Correlation with solar activity 

3.2.1. Solar cycle variation 

The temporal variations in frequency shifts with various measures of the solar activity 
are plotted in Figures 4 and 5 for GONG and MDI data sets, respectively. Although the 
frequency shifts follow the general trend of the solar activity, there are variations between 
different activity indices. For example, we note discrepancies between Su, and Rj, MWSI 
and TSI at high activity period. We examine this variation by calculating the correlation 
between them. This is displayed in Table 3, where linear (rp) and rank (rg) correlation 
coefficients between frequency shifts and the measures of solar activity are tabulated. The 
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linear correlation coefficient is an estimate of the linear relationship between the two sets 
and the rank correlation shows whether one set of numbers has any effect on another set of 
numbers. The high values of rp (> 0.9) clearly indicate a linear relationship between two 
quantities. Therefore, we fit the shifts against activity using linear expression 

5v = a + bl (3) 

and investigate the sensitivity of frequency shifts on solar activity. The slope, b, which 
measures the shift per activity index, and intercept, a, obtained from the linear least-square 
fit for each activity index, J, are also given in Table 3. The gradients obtained from the linear 
regression analysis show a higher sensitivity in the GONG data as compared to the MDI 
and is consistent with the notion that the uncertainty in the determination of frequencies is 
inversely proportional to the square root of the length of the time series. 

A comparison of correlation coefficients for all activity indices indicates that the best 
correlation is obtained for chromospheric activity indices, in particular F w . The recent 
analysis of low-degree modes also report better correlations fo r F 10 , He I equivalen t width and 



Mg II core-to- wing ratio for three consecutive cycles 21-23 (jChaplin et al.l 120071 ). However, 
these authors report relatively low correlation coefficients for these indices as compared to 
our analysis, although the length of the time samples in both cases is the same, i.e. 108 
days. The low-degree modes travel deep to the solar interior while the intermediate degree 
modes are confined to the radiation and convection zones. Virtually all of the activity effects 
on the modes is at the surface. The surface spatial patterns of the low-£ modes are not well 
matched to the spatial distribution of surface activity, hence one may expect intermediate- 
and high-degree modes to correlate better with the activity. We also obtain higher correlation 
co efficients for the 108- day GONG data sets comp ared to the 72-day MDI data. The analysis 



of I Chaplin et al.l (120011 ) and lTripathy et al.l (120071 ) on many different time scales show similar 



correlation patterns. 

Although cycle 23 has been a weak cycle compared to its predecessor, it recorded two 
prominent peaks during solar maximum. The first peak in sunspot number (mid-2000) was 
higher than the secondary peak (late-2001), while several other solar indices (e.g., KPMI, 
Fio, Mg II) had a maximum at the secondary peak. The frequency shifts also exhibited the 
maximum during the secondary peak. Thus R\ is the only activity proxy whose maximum 
does not coincide with the maximum of the frequency shifts. 
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3.2.2. Phase-wise variation 

In order to study the correlation between frequency shifts and activity at different phases 
of the solar cycle, we divide the activity cycle into three phases; the rising activity period 
from 1996 September 22 to 1999 June 26 {Phase i), the high-activity period from 1999 June 
27 to 2003 January 12 (Phase II), and the declining activity period from 2003 January 13 to 
2007 July 26 (Phase III). It should be noted that the relatively flat portion of solar activity 
corresponding to the beginning of minimum phase is not included in this analysis. The start 
and end dates in each activity phase are the same for both GONG and MDI data sets. The 
number of data sets representing Phase I, II and III are; 26 (GM 16 - 41), 34 (GM 44 - 77) 
and 44 (GM 80 - 123) for GONG and 11, 17 and 22 for MDI. Our analysis is restricted to 
Phases / and //for indices KPMI and He I due to the unavailability of activity data during 
the entire declining phase. 

The calculated linear correlation coefficients, rp, for both GONG and MDI data sets 
for all the three phases are shown in Figure 6. In the same figure, we also compare these 
coefficients with those obtained for the complete cycle. It is interesting to note that the 
correlation between frequency shifts and solar activity changes significantly from phase to 
phase except for Fi . In all cases, the rising and the declining phases are better correlated 
than the high-activity phase, except for CI where we obtain a decrease in correlation during 
the declining phase for the frequency shifts calculated from GONG. A substantial decrease 
in correlation during high-activity phase is obtained for R\, MWSI, TSI and He I indices. 
This poor correlation is also visible in Figures 4 and 5 where variation in Sis shows a different 
pattern than the activity proxies. We point out that the R\, MWSI and TSI indices are most 
affected by the sunspots with strong magnetic fields, while other proxies (e.g. MPSI, Mg II) 
are sensitive to the weak component of the magnetic field (see Table 2). The KPMI and 
F w indices are modulated by both components, however the dominant contribution to F w 
comes from the weak component. The change in correlation coefficients for He I is similar 
to that we obtain for CI, but the coefficients for the declining phase can not be compared 
due to insufficient measurements of the He I index. It should be noted here that these two 
indices are influenced by the changes occurring above the chromosphere. 

Because the correlation coefficients clearly reveal a phase-wise variation, we further in- 
vestigate their relation to the corresponding linear regression slopes. The slopes for different 
phases calculated using Equation 3 are tabulated in Table 4. It is seen that the slopes in the 
rising and falling phases are greater than that obtained for the high activity period. How- 
ever, KPMI appears to be an exceptional case where the slopes in the high activity period 
are marginally higher (within la error) than the rising activity period, for both GONG and 
MDI data sets. 
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The calculated linear correlation coefficients between phase-wise correlation coefficients 
and slopes (rp([rp(I,5u)],[b])), and solar activity (rp([rp(I,Su)],[I])) for different activity 
indicators are summarized in Table 5 for both GONG and MDI data sets. In all cases, 
except for CI, we obtain a significant positive correlation between correlation coefficients 
and slopes, and a negative correlation between correlation coefficients and activity indices. 
In general, we find that the GONG data has a hig her correlation t h an th e MDI data. Our 
results are further consistent with the findings of iTripathy et al.l (120071 ) where year-wise 
variations in the linear-regression slopes and the correlation coefficients are measured from 
frequencies that were determined from time series of 9 days. They also found a positive 
correlation between these two but no significant correlation or anti-correlation between the 
slopes and the solar activity. This difference may be attributed to the fact that the frequency 
shifts track the changes in solar activity much better on longer time scales, as discussed in 
the previous subsection. 



3.3. Comparison between rising and declining phases: the hysteresis patterns 



As shown in Section 3.2.2, there is no significant difference between correlation co- 
efficients in rising and declining phases, except for CI, one may expect that the curves 
of frequency shifts plotted against solar activity during the two separate phases should 
overlap with each other or have a constant shift between them. However, the magnetic 
field indices in cycle 22 have been found to foll ow different curves during ascending and 



desce nding phases, exhibiting a hysteresis loop (jJimenez-Reyes et al.l Il998t ITripathy et al. 



2001 ). It is also shown tha t not all activity indices have a linear relationship among them 
( Bachmann fc White I1994J ). Since cycle 23 has been referred to as anomalous solar cycle 
( jde Toma et al.l 12004 ) due to differences in various activity measures as compared to other 
cycles, we plot frequency shifts against the activity and examine the behavior of these curves 
to obtain a deeper insight of the variation of frequencies with solar activity. 

In Figures 7 and 8, we show the variation of 8v as a function of four different activity 
indices (F 10 , KPMI, MPSI and Mg II) for GONG and MDI data sets, respectively. The 
open and filled symbols represent the ascending and descending phase, respectively; the 
high activity period is shown by asterisk symbols. It is evident that the mode frequency 
shifts follow similar paths during the ascending and descending phases for the radio flux, 
F w , and full-disk magnetic field, KPMI. However, the other two pairs (MPSI, and Mg II) 
exhibit different paths and form a broader curve. The behavior of Fiq, MPSI and Mg II in 
cycle 23 is consistent with that for cycle 22 while the vari ation of frequencies w ith KPMI 
contradicts the findings of I Jimenez- Reyes et al.l (119981 ) and ITripathy et al.l (120011 ). We note 
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that the instrument at Kitt Peak was upgraded in 1992 to reduce the noise level. The post- 
upgrade data were calibrated to place them at the same scale as pre-upgrade data but any 
systematic offset between the two data sets could lead to a constant shift in the values of 
KPMI and may account for the hysteresis pattern seen in cycle 22. 

4. Discussion 

The intermediate-degree modes are generally confined to the convection zone, therefore 
one may expect the variation in their characteristics to be influenced by changes in this 
region. The solar cycle related shifts are believed to occur primarily at the surface reflection 
of the modes. In addition, their variation and strong correlation with solar activity indices 
over the solar cycle provide sufficient evidence for their association with the changing solar 
magnetic activity. Various measurements have been made to estimate such changes in mag- 
netic activity. In this analysis, we have considered three proxies which represent quantitative 
measurements of the magnetic field strength, i.e. KPMI, MWSI and MPSI. As outlined ear- 
lier, the KPMI is total line-of-sight magnetic field while MWSI provides only the strong 
component in the photosphere. On the other hand, MPSI estimates weak fields associated 
with magnetic plages in the chromosphere. 

In order to study the differences in the correlation between frequency shifts and various 
magnetic measurements, we plot 108-day averaged photospheric fields against the MPSI 
(Figure 9). Note that the weak and strong fields do not vary in a similar way at all activity 
phases. We find a linear relationship during rising and falling activity periods among all three 
proxies, but a large deviation from the linear trend is seen between MPSI and MWSI during 
the high-activity period. The variation in the KPMI is linear with the MPSI in rising- and 
high-activity periods, however, their behavior remains unexplored during the declining phase 
due to the unavailability of KPMI data. To supplement the KPMI measurements during the 
declining phase, we include in the same figure the line-of sight magnetic field measured by 
Synoptic Optical Long-term Investigations of the Sun (SOLIS) instrument located at Kitt 
Peak. The SOLIS observation^ are made using a different spectral line (Fe I 630.2 nm) than 
the KPMI, but the overlapping data of these two instruments show a similar response to 
active regions and other solar phenomena (J. W. Harvey, private communication). However, 
we notice a downward shift in the SOLIS fields as compared to the KPMI which may be due 
to the use of a different spectral line. Moreover, we find a strong correlation between SOLIS 
magnetic field and the MWSI which is consistent with the results obtained for KPMI. 



3 http://solis. nso.edu/vsm/svsm_ml l_meanfield.dat 
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Therefore, we do not expect frequency shifts to correlate equally with all activity indi- 
cators at different phases. Our analysis for the complete cycle and different activity phases 
clearly shows that the frequency shifts have the best correlation with the 10.7 cm radio flux. 
The radio flux is measured in the high chromosphere and lower corona, but it is modified 
by the evolution of sunspots (strong magnetic field) and radio plage s (weak magnetic field). 
By analyzing low-£ frequencies for cycle 21-23, 1 Chaplin et al.l (120071 ) concluded that the fre- 
quencies have good sensi t ivity to the effects of the weak component of the magnetic fields. 
Recently, iTripathy et ajj (120071 ) have shown that activity related changes in oscillation fre- 
quencies can be tracked on a time scale as short as nine days. They also obtained the best 
correlation with the radio flux. In a theoretical study, ICampbell &: Roberts I (119891 ) also 
showed that changes in the chromospheric magnetic field strength can cause changes in the 
2?-mode frequenc ies at higher I { > 100 ). On the other hand, using the technique of local 
helioseismology, iHindman et al.l (120001 ) argued that the magnitude of the variation in fre- 
quencies depend upon the magnetic field strength associated with active regions on the solar 
surface and the physical phenomenon inducing the shifts is confined to the surface layers 
of the Sun. On the basis of our study, we find that both components of the magnetic field 
influence the changes in the oscillation frequencies, however the weak component plays a 
dominant role. To examine this conclusion, we developed various regression models between 
MWSI (strong magnetic field), MPSI (weak magnetic field), and F\§ and bv from GONG 
and MDI data sets. The variance for different regression models calculated using frequency 
shifts from GONG data are given Table 6. It is clearly seen that we obtain higher correlation 
at all activity phases for the models based on either F w or combined MPSI and MWSI data. 
These models are able to explain more than 90% of the variation in the frequencies at all 
activity phases. In contrast, the model with strong component of magnetic field (MWSI) 
explains less than 50% of the variation while the model with only weak component (MPSI) 
explains less than 80% of the variation in frequencies at high-activity period. A similar 
behavior is also seen in case of MDI data. 



5. Summary 

The improved and continuous measurements of intermediate-degree p-mode frequencies 
for a complete solar cycle demonstrate that, while the frequencies vary in phase with the 
solar activity, the degree of correlation between the frequency shifts and activity indices 
differs from one phase of activity to another. 

Although there is a strong correlation between frequency shifts and all activity proxies 
during the rising and declining phases, we find a significant change in the correlation at 
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the high-activity period. The maximum decrease is obtained for Ri, MWSI, TSI, and He I, 
where first three indices are most affected by the sunspots with strong magnetic fields in the 
photosphere while He I is influenced by the change in temperature in transition region. The 
change in correlation coefficients for CI is similar to that of He I but we could not compare 
the coefficients in the declining phase due to insufficient measurements of the He I index. We 
further find a negligible change in the correlation for the 10.7 cm radio flux which provides 
estimates of the variation in both weak and strong components of the magnetic field. It 
is also influenced by a slowly varying component (S'-component) that originates thermally 
in localized regions of high electron densities and magnetic fields present in the vicinity of 
sunspot and chromospheric plages. These regions remain present over several rotations of 
the Sun and contribute to the change in oscillation frequencies. Our analysis suggests that 
both components of the magnetic field contribute significantly to the variation in oscillation 
frequencies, however the contribution from the weak component appears to be dominant. 

The correlation coefficients are strongly correlated with the slo pes obtained fr o m the 



linear regression analysis, in agreement with the earlier findings of iTripathy et al.l (120071 ) 



using frequencies determined from time-series of 9 days in length. In addition, slopes and 
correlation coefficients are found to be strongly anti-correlated with the measures of solar 
activity. However, the coronal index shows a poor correlation in all cases. It may be due 
to the fact that the solar oscillations are most affected by the magnetic activity of lower 
atmosphere and are not influenced by the changes in corona. 
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Fig. 1. — Duty cycles of (a) GONG and (b) MDI time series. 




Fig. 2. — Comparison of common p modes observed in all (a) GONG, (b) MDI and (c) 
combined GONG and MDI data sets. 
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Fig. 3. — Temporal evolution of frequency shifts calculated from the GONG (open symbols) 
and MDI (filled symbols) frequencies. The dashed line represents the scaled 10.7 cm radio 
flux. 
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Fig. 4. — Temporal evolution of the GONG frequency shifts (symbols) with different activity 
proxies (filled regions). 
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Fig. 5. — Same as Figure 4 but for the shifts calculated from the MDI frequencies. 



-20- 




Fig. 6. — Bar-chart showing phase-wise variation of the Pearson's linear correlation coeffi- 
cients for different activity indices. The complete solar cycle is divided into three phases; 
the rising (I), high (II) and declining (III) activity Each activity phase has values for the 
GONG (left) and the MDI (right). Note that missing values for KPMI and He I are due to 
unavailability of activity measurements during the declining phase. 
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Fig. 7. — Variation of GONG frequency shifts with different activity proxies. The rising, 
high and falling phases are shown by open symbols, asterisks and filled symbols, respectively. 
The symbols have been joined by solid line to guide the eye. 
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Fig. 8. — Same as Figure 7 but for the MDI frequency shifts. 
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Fig. 9. — (left) Variation in KPMI (total field) and MWSI (strong field) as a function of 
MPSI (weak field). Rising, high and falling activity periods are shown by the squares, filled 
circles and pluses , respectively. The crosses are for the SOLIS data. The solid lines denote 
the linear fits at different phases, (right) The linear correlation coefficients, rp (filled) and 
rank correlation coefficients, (hatched) at different phases corresponding to the left panels. 
In the lower panel, rp and rs between the SOLIS measurements and MPSI are shown by 
open and diagonally hatched bars, respectively. 
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Table 1. Description of frequency data sets. 





GONG 


MDI 


Period 


1996 May 1 - 2007 Aug 31 


1996 May 1 - 2008 Feb 27 


Number of datasets 


113 


58 


Extent (days) 


108 


72 


Frequency range (mHz) 


1.5 < v < 4.0 


1.5 < v < 4.0 


Number of common modes 


500 


885 


Degree range 


19< I <147 


7< I <188 


Radial nodes 


1< n <19 


1< n <21 









Table 2. Various proxies for solar activity. 



Activity 


Spectral line 


Influential 


End date a 


Data source 


Index 




field component 






Ri 




Strong 


2008 27 Feb 


|ftp://ftp.ngdc.noaa.gov/STP/SOLAR.DATA/SUNSPOT.NUMBERS/| 


KPMI 


Fe I 868.8 nm 


Both 


2003 Sep 21 


ftp: / / nsokp.nso.edu /kpvt /daily/stats / mag.dat 


MWSI 


Fe I 525.0 nm 


Strong 


2008 27 Feb 


ftp:/ /howard. astro, ucla.edu/pub/obs/mpsi_data/index.dat 


Fio 




Both 


2008 Feb 27 


|ftp://ftp.ngdc.noaa.gov/STP/SOLAR.DATA/SOLARJlADIO /FLUX/| 


He I 


He I 1083.0 nm 




2003 Sep 21 


ftp : / /nsokp . nso . edu /kpvt /daily/stats/he.dat 


MPSI 


Na I 589.5 nm 


Weak 


2008 27 Feb 


ftp:/ /howard. astro, ucla.edu/pub/obs/mpsi_data/index.dat 


Mg II 


Mg II 279.9 nm 


Weak 


2007 Oct 6 


http: / /www. sec. noaa.gov/ftpdir/sbuv/NOAAMgll.dat 


TSI 






2008 Feb 27 


ftp:/ /ftp.pmodwrc.ch/pub/data/irradiance/virgo/TSI 


CI 


Fe XIV 530.3 nm 




2005 Dec 31 


|ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_CORONA/| 



a Start date for all data sets is 1996 May 1. 
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Table 3. The correlation coefficients and the linear regression parameters between activity 
indices (/) and frequency shift (Su) for the total period covered in this analysis. Shown 
here are the Pearson's linear coefficient (rp), Spearman's rank correlation (rs), the 
intercept (a) and the mean shift per unit change in activity (6). The two-sided significance 
(Ps) is smaller than 10~ 16 and is not shown here. 



Activity Index 


Data Source 


r P (I,Su) 


r s (I,5v) 


a 

( M Hz) 




b 

(/iHz per activity unit) 


Ri 


GONG 


0.98 


0.98 


-6.28X10- 2 


± 


1.60x10" 


-4 


9.07X10- 4 ± 2.28xl0- 6 




MDI 


0.97 


0.95 


-4.01X10- 2 


± 


2.25x10" 


-4 


8.31X10- 4 ± 3.22X10- 6 


KPMI 


GONG 


0.99 


0.96 


-1.11x10-2 


± 


4.09x10" 


-4 


8.63X10- 3 ± 2.87X10- 5 




MDI 


0.98 


0.93 


-8.47x10-2 


± 


5.92x10" 


-4 


8.02X10- 3 ± 4.16X10- 5 


MWSI 


GONG 


0.95 


0.95 


-5.54x10-2 


± 


1.48x10" 


-4 


1.83X10- 1 ± 4.75X10- 4 




MDI 


0.95 


0.94 


-3.43x10-2 


± 


2.10x10" 


-4 


1.67X10- 1 ± 6.59X10- 4 


Fio 


GONG 


0.99 


0.98 


-1.17X10- 1 


± 


2.80x10" 


-4 


8.76X10- 4 ± 2.18X10- 6 




MDI 


0.98 


0.97 


-9.11x10-2 


± 


3.95x10" 


-4 


8.14X10- 4 ± 3.10X10- 6 


He I 


GONG 


0.97 


0.95 


-1.84X10- 1 


± 


6.59x10" 


-4 


2.89X10- 3 ± 9.81 xlO- 6 




MDI 


0.96 


0.95 


-1.55X10- 1 


± 


9.58x10" 


-4 


2.71X10- 3 ± 1.42 xlO- 5 


MPSI 


GONG 


0.98 


0.97 


-5.98x10-2 


± 


1.54x10" 


-4 


3.97x10-2 ± l.OOxlO- 4 




MDI 


0.98 


0.97 


-3.83x10-2 


± 


2.19x10" 


-4 


3.71x10-2 ± 1.43xl0- 4 


Mg II 


GONG 


0.95 


0.94 


-1.92 


± 


4.96x10" 


-3 


7.01 ± 1.82 xlO" 2 




MDI 


0.95 


0.93 


-1.75 


± 


7.22x10" 


-3 


6.49 ± 2.65 XlO" 2 


TSI 


GONG 


0.91 


0.89 


-108.83 


± 


2.95x10" 


-1 


7.97x10-2 ± 2.16xl0- 4 




MDI 


0.91 


0.90 


-104.31 


± 


4.32x10" 


-1 


7.64x10-2 ± 3.16X10- 4 


CI 


GONG 


0.96 


0.94 


-8.37x10-2 


± 


2.69x10" 


-4 


1.18x10-2 ± 3.59xl0- 5 




MDI 


0.94 


0.92 


-5.84x10-2 


± 


3.88x10" 


-4 


1.08x10-2 ± 5.17xl0- 5 
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Table 4. Calculated linear regression slopes (in /xHz per activity unit) for different activity 
phases of solar cycle. Note that the slopes for the total period are different from Table 3 as 
these are calculated for a different period. Here periods covered for GONG and MDI are 
same in all activity periods. Missing values of KPMI and He I for Phase III are due to the 

unavailability of activity measurements. 



Activity 
Index 


Data 
Source 


Phase I 




Phase II 




Phase III 






Total 




Ri 


GONG 


8.56X1CT 4 


± 


4.27x10" 




5.72xl0- 4 


± 


1.12x10" 


-4 


9.26x10 


" 4 ± 3.74x10- 


-5 


9.07x10- 


-4 


± 1.74x10" 


-5 




MDI 


6.66X10- 4 


± 


7.39x10" 


-5 


5.56X10- 4 


± 


1.62x10" 


-4 


8.82x10 


" 4 ± 9.43x10- 


-5 


8.35x10- 


-4 


± 3.17x10" 


-5 


KPMI 


GONG 


8.56X1CT 3 


± 


2.72x10" 


-4 


8.69xl0- 3 


± 


9.49x10" 


-4 








8.77x10- 


-3 


± 1.66x10" 


-4 




MDI 


7.40X10- 3 


± 


8.19x10" 


-4 


7.96X10" 3 


± 


1.68x10" 


-4 








8.06x10- 


-3 


± 3.35x10" 


-4 


MWSI 


GONG 


2.63xlCr 1 


± 


1.33x10" 


-2 


7.65X10- 2 


± 


1.57x10" 


-2 


1.68x10 


- 1 ± 7.99x10- 


-3 


1.79x10" 


-1 


± 5.86x10" 


-3 




MDI 


2.19X10" 1 


± 


1.95x10" 


-2 


8.28x10-2 


± 


2.14x10" 


-2 


1.59x10 


- 1 ± 1.98x10" 


-2 


1.66x10" 


-1 


± 7.97x10" 


-3 


Fio 


GONG 


9.34X1CT 4 


± 


3.10x10" 


-5 


6.27xl0- 4 


± 


3.24x10" 


-5 


9.54x10 


" 4 ± 2.34x10" 


-5 


8.71x10" 


-4 


± 1.12x10" 


-5 




MDI 


7.82X1CT 4 


± 


6.32x10" 


-5 


6.16xl0- 4 


± 


7.71x10" 


-5 


9.51x10 


" 4 ± 7.29x10" 


-5 


8.10x10" 


-4 


± 2.16x10" 


-5 


He I 


GONG 


2.99X1CT 3 


± 


1.57x10" 


-4 


1.52xl0~ 3 


± 


2.96x10" 


-4 








2.88x10" 


-3 


± 9.80x10" 


-5 




MDI 


2.57xl(r 3 


± 


2.89x10" 


-4 


1.58X10- 3 


± 


4.71x10" 


-4 








2.69x10" 


-3 


± 1.45x10" 


-4 


MPSI 


GONG 


4.61xl(T 2 


± 


1.42x10" 


-3 


2.44x10-2 


± 


2.65x10" 


-3 


4.16x10 


"2 ± 9.54x10" 


-4 


3.96x10" 


-2 


± 7.92x10" 


-4 




MDI 


4.02X10- 2 


± 


3.40x10" 


-3 


2.31x10-2 


± 


4.48x10- 


-3 


4.15x10 


"2 ± 2.32x10" 


-3 


3.69x10" 


-2 


± 1.21x10" 


-3 


Mg II 


GONG 


7.14 ±2.15x10" 


-1 


4.31 ±4.98x10" 


-1 


8 


.03 ±2.06x10- 


-1 


7.10 ±2.15x10" 


-1 




MDI 


6.08 ±5.13x10" 


-1 


4.06 ±9.23x10- 


-1 


7.93 ±5.25x10" 


-1 


6.60 ±3.14x10" 


-1 


TSI 


GONG 


9.47X10- 2 


± 6.12x10" 


-3 


1.50x10-2 


± 7.60x10- 


-3 


7.48x10" 


"2 ± 6.81x10" 


-3 


8.03x10" 


-2 


± 3.74x10" 


-3 




MDI 


8.03x10-2 


± 1.26x10" 


-2 


1.56x10-2 


± 1.18x10" 


-2 


6.94x10 


" 2 ± 8.91x10" 


-3 


7.52x10" 


-2 


± 5.26x10" 


-3 


CI 


GONG 


9.46X10- 3 


± 4.78x10" 


-4 


7.87xl0- 3 


± 1.10x10" 


-4 


9.33x10 


~ 3 ± 2.47x10- 


-3 


1.21x10" 


-2 


± 3.92x10" 


-4 




MDI 


7.95 xlO- 3 


± 


8.44x10" 


-4 


6.65xlQ- 3 


± 1.97x10" 


-4 


1.04x10- 


" 2 ± 3.01x10- 


-3 


1.11x10" 


-2 


± 6.18x10" 


-4 
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Table 5. The correlation statistics between phase-wise correlation coefficients and slopes 
(rp([rp(I,5u)},[b])), and solar activity (rp([rp(I, Su)],[I])) for different activity indicators 

using GONG and MDI data sets. 



Activity Index 


rp([r P (I,8»Mb]) 


rp([rp(I,8u)],[I]) 


GONG 


MDI 


GONG 


MDI 


Ri 


0.98 


0.74 


-0.92 


-0.85 


KPMI 


0.90 


0.74 


-0.85 


-0.77 


MWSI 


0.85 


0.94 


-0.95 


-0.85 


fio 


0.93 


0.87 


-0.89 


-0.72 


He I 


0.99 


0.99 


-0.85 


-0.82 


MPSI 


0.97 


0.97 


-0.96 


-0.92 


Mg II 


0.87 


0.97 


-0.98 


-0.95 


TSI 


0.99 


0.99 


-0.85 


-0.87 


CI 


0.47 


0.34 


-0.05 


-0.53 



Table 6. Calculated variance for different regression models using GONG frequency shifts 
(Su) and activity indices (Fi , MPSI and MWSI) for different phases of the solar cycle. 
The coefficients A, B and C are determined from the multiple linear regression fit. 





Regression Model 


Total 


Phase I 


Phase II 


Phase III 


8v 


= A- 


h BF 10 


0.99 


0.98 


0.96 


0.99 


5u 


= A- 


- B (MPSI+MWSI) 


0.97 


0.99 


0.92 


0.99 


8u 


= A- 


- B MPSI+ C MWSI 


0.97 


0.99 


0.93 


0.99 


8u 


= A- 


r B MPSI 


0.96 


0.99 


0.77 


0.98 


8v 


= A- 


h B MWSI 


0.90 


0.93 


0.47 


0.84 



